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In this paper we show that (a) all the known exact solutions of the problem of N- 
anyons in oscillator potential precisely arise from the collective degrees of freedom, 
(b) the system is pseudo-integrable ala Richens and Berry. We conclude that the 
exact solutions are trivial thermodynamically as well as dynamically. 
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I. INTRODUCTION 



Non-relativistic quantum mechanics in two space dimensions admits the possibility of 
fractional statistics and particles obeying fractional statistics are known as "anyons". 
Anyons are objects defined in the quantum framework described by a Lagrangian of the 
form, 

L = \ E rl + a £ % - V{?& 6 l3 = tan^f^-, (1) 

z i=i i<j Xl x i 

where r*j denote particle coordinates and V(fi) is some confining potential which we choose 
to be the harmonic oscillator potential V(fi) = J2iL± ■ The a dependent term is the 
statistical interaction. The harmonic oscillator potential is convenient since the dynamics 
of the system is well understood in the limit a = . We may choose to describe the 
system of anyons either through multivalued wavefunctions which occur naturally in the 
quantum mechanics on multiply connected spaces or equivalently interms of single valued 
wave functions in the presence of statistical interaction. Quantum mechanically a system 
of N-anyons confined in a harmonic oscillator potential is an exceptional system since many 
exact solutions to the energy eigenvalue problem are known even when the many body 
Hamiltonian is nonseparable. 

The known results about the spectrum of N-anyons confined in oscillator potential fall 
into two categories: (a) Exact eigenvalues which are linear functions of the statistical pa- 
rameter a, Ej m = 2m + \j — a N \ N ^~ 1 > | + N, where j and m denote the angular and radial 
quantum numbers, and (b) eigenvalues that are nonlinear functions of a as evidenced by 
numerical calculations for iV = 3, 4 anyons p|-4T0|] and also meanfield calculations for large iV 



||1 1|| . The special case in this context is iV = 2 where the spectrum is completely subsumed 
by the category (a). The nonlinear spectrum displays many level crossings, some of which 
are of the Landau-Zener type(avoided crossings). [jSHICfl This has lead to the conjecture that 
the anyon system may be nonintegrable Jl0| ( or even chaotic) for N > 3. 

Two comments are in order here: Firstly, what is the reason for the existence of exact 
solutions of the type (a) in this nonseparable many body system? Infact the approach of 



Basu etal [[?]] shows that j=0 exact solutions form a subspace of the full Hilbert space which 
block diagonalises the Hamiltonian. This indicates that the system may be atleast partially 
separable. The origin of this has not been elaborated yet to the best of our knowledge. The 
second comment pertains to the "integrability" of the system. If one considers the classical 
Lagrangian, the statistical interaction is a total derivative. The Euler-Lagrange equations 
of motion are therefore the same as that of an 2iV-dimensional oscillator and this ofcourse 
is an integrable system. Why do then the numerical results suggest nonintegrability? 

In this paper we elaborate on and answer the above questions. In section 2, we briefly 
mention the known exact solutions to the quantum many body problem. In section 3 we 
outline the classical Lagrangian formulation of N-anyons and prove the system is "integrable" 
in the Liouville sense, ie., it admits 2N constants of motion in involution. In section 4 we 
exhibit a choice of coordinates at classical level which gives partial separability. From this we 
can isolate a set of collective coordinates whose semiclassical quantisation leads to the known 
exact solutions. These solutions correspond to the periodic orbits in the phase space of the 
oscillator which continue to remain periodic even in the presence of anyonic interaction. In 
section 5 we consider the integrability aspect of the system. We argue that even though we 
have 2N constants of motion in involution there exist invariant surfaces which do not have 
the topology of 2N dimensional torus. Following Berry-Richens [13| we conclude that the 
many anyon system confined in oscillator potential is pseudointegrable. 
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II. EXACT SOLUTIONS 



The quantum mechanical Hamiltonian corresponding to the Lagrangian in Eq.(l) is given 

by, 

1 N 1 N N £■■ a 2 N r- r-u 

H=ilEti + - 2 Er>-a E 3 + (2) 

^ z=l ^ i=i j>i=i 1 ij * ijtj,k'ij'ik 

where 

4j = - ij-) x (Pi -ft-)- 

and all distances have been expressed in units of 1/y/mu), where m is the mass of the particle 
and u) is the oscillator frequency. Notice that the statistical interaction is independent of the 
centre of mass. We now briefly discuss the class of exact solutions already known. 0-0] For 
discussing the known class of exact solutions it is convenient to use the complex coordinates 
Zi, z~i in terms of which the Hamiltonian takes the form, 

# = -2£^4E^-«E(^-^) + y E-^ 5 (3) 

i ^ i i<j Z V z ijtj,k ^j z ^k 

where di = d/dzi] dij = di — dj, etc., and the eigenvalue equation is, 

Hip(zi, = Eip(zi, z^. (4) 
The conserved angular momentum J with eigenvalues denoted by j , is given by, 

J = Zidi - Zidi. (5) 

Defining variables, 

N N 

x = H z 'iji t = ^2ziZi, (6) 

i<j i=l 

we may classify the known exact solutions as follows: 
(a)j < solutions- 

^ = ixr^e-*/ 2 ; j < (7) 



with the energy eigenvalues given by, 

E=N - j+a mzA. (8) 

Here and in what follows <f>j generically denotes an eigenfunction of the angular momentum 
operator J with eigenvalue j. For elucidating the solutions the specific form of is irrelevent. 

(b) j = solutions- 

= \X\ a Mt)e- t/2 ; J = 0. (9) 

with (j) (t)(= J2T=i Ckt k ) a polynomial of degree m in t. The corresponding energy eigenval- 
ues are given by, 

E = N + 2m + a N ( N ~ 1 \ (10) 

The second solution is necessarily bosonic since t is symmetric, where as the first solution 
needs explicit symmetrization and antisymmetrization of the wavefunction in terms of 2i to 
obtain the bosonic and fermionic wavef unctions. Since this is always possible, the degeneracy 
of the first type solution is exactly the same for both bosonic and fermionic type solutions 
for any given angular momentum j (< 0). We may also take a combination (product) of 
the solutions the types discussed above, to get further j < solutions. This then generates 
a infinite tower radial excitations for each value of j. 

(c) j > solutions- 

^ = |Xr<^)e-'/ 2 ; j > (11) 
with the energy eigenvalues given by, 

E = N+J . a !^. (12) 

Caution must be excercised in choosing the value of j for these solutions since the wave 
function is not square integrable for all values of j. Infact the lower bound is obtained 
requiring that the wave function be square integrable over the whole domain of a(0 < a < 1). 



This means that j > (N — 1)(N — 2)/2. If however this condition is not satisfied then the 
wave function remains regular only for some values of a (0 < a < 2j/N(N — 1)) but not for 
all < a < 1 which gives rise to the so called noninterpolating solutions which have also 



been discussed in the literature. 12 



All these solutions for the energy eigenvalues have a linear dependence on a with a 
coefficient ±iV(iV — l)/2 while the corresponding eigenfunctions are finite order polynomials 
apart from the overall | X \ ±a and the Gaussian factors. These solutions (a)-(c) cover all 
the known exact solutions. However it is by now known that these exact solutions form 
only a subset of the full Hilbert space and existence of nonlinear solutions has been shown 
numerically as well as through meanfield calculations. It is our aim here to understand the 
reason for the existence of this dichotomy. 

III. CONSTANTS OF MOTION 

We begin with the analysis of the classical Lagrangian given by Eq.(l). It is convenient 
to write the Lagrangian in the form, 

L = lj2f,-r1]+aE r -f^p 1 ^ (13) 

z 8=1 i<j 

where the dots indicate the time derivatives. The first step is to introduce relative coordi- 
nates and separate the trivial centre of mass degree of freedom. To this end we write, 



1 a a 

Pa = [ f- = E ^ - \\ -TTVili a = l,...,A-l (14) 

la(a + l)fe=i V a + 1 



with inverse relation given by 



— + E 1 ^ ] + Rem = A1p a + R cm , (15) 



k =i ^k(k + iy 

where p a ,a = 1,...,N — 1 are dimensionless relative coordinates and R cm is the centre of 
mass coordinate. It follows that, 



It is straight forward to see that, 

L = Lqm + L re i, 

where 



1 ^ 2 

Lcm = -j[ R CM ~ r cm}' 



t * r -h. -ai i A ij A ijP a x p 5 

2 i<j ^ij^ijPc-Pd 

where Afj = Af — Aj. From now on we concentrate only on the L rei and drop the subscript. 
It is easy to see that the Euler-Lagrange equations of motion are, 

Pa = -Pa- (16) 

There is no a dependence in these equations since the corresponding term is a total deriva- 
tive. From these equations it follows that 

£a=\\p a + fi]; la = Pa><'p a ; a=l,...,N-l (17) 

are constants of motion. In the Hamiltonian formulation the p a are expressed in terms of 
conjugate momenta and coordinates which do contain the a dependence, 

P - — = p -aA ■ A = V A % A %P*v (18) 

1 ax r% . Hax "^OX! ^ax / y Ar a j -, -, ■> V u / 

dpax X£j A ij A ijPc-Pd 

P — - ' U A • A - \ " A ij A ijP bx /-|q\ 

ay ~ fin ~ P ay a,A -ay, Slay ~ Ac Ad J => ' V iy / 

OPay iKj • kjPc-Pd 

and the constants of motion in relative coordinates are, 

£a = \[{Pax + uAaxf + (Pay ~ OiA ay f + f a \ (20) 
l a = p a X P a - a(p ax A ay + payA ax ). (21) 

Clearly the Hamiltonian is given by 

JV-1 
a=l 
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Note that the a-dependent term in L is well defined only if ^ for all i,j. That is only if 

A%A%p a .p h ± V i,j 

Consequently the expressions for £ , l a and H are also valid only if is not zero. In effect 
the classical configuration space on which L is well defined is the space 

Q N -i = _ A; A = { p a I A c l3 A% p c .p d = for some patr(s) (22) 

The space Q N ~ l is not simply connected. Its fundamental group is the same as the funda- 
mental group of {R 2 — (N — 1) points} which is known to be nontrivial. For N = 2, 7Ti = Z 
while for N > 3, 7Ti is nonabelian. 

The corresponding phase space is the cotangent bundle of Q on which S a , l a and H are 
well defined. This phase space is topologically nontrivial and for our purposes we do not 
need the full machinery for handling this topologically nontrivial phase space. It is sufficient 
to pretend that the full space is still R^' 1 but simply avoid coincident points. 

It is straight forward then to prove the following Poisson bracket relations, 

{£ a ,£ b } = {l a ,l b } = {£ a ,l a } = V a,& = l,...,jV-l, (23) 
{H,£ a } = {H,l a } = 0. (24) 

Thus for 4(iV — 1) degrees of freedom P a , p a we have 2(iV — 1) constants of motion in 
involution. So a necessary condition for a system to be integrable is satisfied. We will 
tentatively refer to this system as being integrable (or potentially integrable) in the "Liouville 
sense" . 



IV. COLLECTIVE MODES 

The anyonic interaction term in L is invariant under two sets of time independent trans- 
formations: 

Pa^Va = R(0)p a V a, (25) 

Pa-^rj a = \p a V a, (26) 



where R(9) denotes a rotation of the vector by an angle 9 in two dimensions. Therefore 
for any given configuration {p a } at any given time t we can always rotate the axes so that 
Ply = (say). Thus by making a time dependent rotation we can ensure that i]i y = for all 
t. Therefore define, 

(cosO —sin9\ 
Va, (27) 
sinO cosO J 

The following identities are easy to prove: 

Pa-Pb = Va-Vb, (28) 

p a x p b = rf a x ff b , (29) 
p a = ff a + 26ff a xff a + 2 € (30) 
Pa x p b = 9rf a .rf b + rf a x ff b . (31) 

Now the set {r/ a ; a = 1,...,N — 1} is effectively a 2(iV — 1) — 1 dimensional vector since 
rj 1: y = for all time t. We can therefore introduce the standard spherical coordinates by the 
usual procedure, 



where 



= R£a, £i y — 0; ^2^,1-1, 

a 


(32) 


^N-a,x = SlS2---S2a-2C2a-l, 


(33) 


^N-a,y — SiS2---S2a-lC2a, 


(34) 


?l,x = s l s 2---S2AT-4, 


(35) 




(36) 



where a = 2, N — 1 and s M = sin 6> M , c M = cos 6^. Interms of these variables the Lagrangian 
may be rewritten as 

L = i [A - # + RV- + 2RH Z C x I + R 2 E I] + a^f^H + a £ ^ffi^-fr , 

(37) 



Since fj a are obtained from p a by the same rotation matrix for all a, the angle 9{t) clearly 
describes a collective rotation of all the N-anyons about the centre of mass(say). The anyonic 
interaction term ( a total time derivative) is manifestly independent of R(t). This may also 
be regarded as a collective mode as discussed below. It is the semiclassical quantisation of 
these two modes that yields the exactly known energy eigenvalues. To elaborate these points 
let us consider the Euler-Lagrange equations of motion. Since the a dependent part of the 
Lagrangian is a total time derivative we can ignore it for analysing the classical equations of 
motion. Clearly these equations are identical to that of the oscillator equations of motions 
as given in Eq.(16). Translating this into the spherical coordinates we obtain, 

[R + R- R9% - [R9 + 2R9] vg a + 2[R- R9V]i a + R£ a = 0, (38) 

where 



V 



1 
-1 

and in the matrix equation above the £ is a column vector with two elements. The matrix V 
between two vectors essentially generates their cross product. Taking the dot product with 

— * 

£ a and summing over a we find, 

[R + R- R9*} = 2R9 J2 £ x t + R £ 'C (39) 

a a 

where we have made use of the identities, 

- u - % -J 

a a a 

— * 

. Taking the cross product with £ a and summing over a we find, 

[R0 + 2 R9] = -2R & xt-RE ^ x C (40) 

a a 

Using the above two relations Eq. [38] can be rewritten as, 

iC + 2[R - R9V]l + [2R9% x & + + [2^g x f 6 + x Qv& = 0, (41) 

where sum over b is assumed. Here the equations motion for R and 9 are still coupled to all 
the other internal coordinates and as yet they are not collective coordinates. We therefore 
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need to impose a set of initial conditions which may lead to the separation of these two 
modes as collective modes from the internal variables 9^. To realise this let at time t — 0, 
all velocities £ a (t = 0) =0 for all a. Then the above equations reduce to, 

R[£a + (E&x&)V&]=0, (42) 

b 

R + R - R9 2 = 0, (43) 
R0 + 2R6 = -R^2i b x ( b . (44) 

b 

Now consider the first equation for R ^ and a — 1. From the initial conditions it is 
obvious that, 

&* = 0; 6xE(6 x & = 

because £i y = for all t. 

Now i/£i x is nonzero, then x ^ b = 0, and hence 

6 = V 6. 

Since the equations of motion are second order in t, we have proved that: if ^1(0) = 
V a, and -R(O), £i x (0) are non zero then\/ t, 

R + R- R9 2 = 0, 
R6 + 2i?fl = 0. 

For oscillator (ct = 0) R(0) and £i x (0) being non zero is a special class of initial condition. 
Indeed in general for a second order equation system if both first and second order derivatives 
vanish at any t then this can hold only for some restricted class of coordinate values. For 
the equations of motions considered in the configuration space -R^ -1 the restricted class is 
precisely characterised by R(0) ^ and £i x (0) ^ However when a^0 the configurations 
space is Q N ~ l which does not contain points which have R = and C,\ x = 0. Thus there are 
no restrictions on the initial conditions in Q N ~ l . 
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The above initial conditions amount to freezing the "internal motion" of the anyons. The 
remaining motion is a collective motion described by R(t) and 6{t) which is just a 2(N — 1) 
dimensional oscillator. R(t) being nonzero implies that the angular momentum must be 
nonzero. This explains in what sense R(t) can be interpreted as a collective mode. For 
describing the motion of collective mode the effective Lagrangian is 

L eff = \[R* + RW - i? 2 ] + a^^le. (45) 

This is identical to the Lagrangian in the relative coordinate for a two anyon system with 
a — > a N( * N 2 ~ 1 * ) anc i j n 2(jV — 1) dimensions. Semiclassical quantisation will then reproduce 
all exactly known energy eigenvalues summarised in Sect. 2. In this effective Lagrangian 
picture the only known memory of N resides in the coefficient of a. This can be understood 
by noting that when all the N-anyons are rotated by 2n about the centre of mass they also 
circle each other to pick up the extra phase. 

If on the otherhand we consider R = 9 = 0aXt = class of initial conditions then we see 
that R, 9 are nonzero and hence R(t), 6{t) do depend on 6^. Thus the collective motion is not 
fully decoupled from the "internal motion". We therefore refer to this system as partially 
separable. 

Incidentally the same conclusions can be drawn from the Hamiltonian formulation. For 
completeness we give the relevent expressions and arguements. The conjugate momenta for 
the hyperspherical variables are given by, 

Pr = R, (46) 
Pe = RH + R 2 Y, 0,F, + a N(N ~ ^ , (47) 
P M = i? 2 rj^ + R 2 QF ii + aG,, (48) 

where, 

r» = sis 2 --s (U _i, fi = 1, ...,2N - 3, n = 1 
and F M and are defined through, 

2N-A N-l 

E w = E £« x H 

fi=l a=l 
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2AT-4 Aa Ab t v t 

2. "mm = lv , c ^ 7 7 • 
The Hamiltonian is then given by 

H = H 1 + H 2 , (49) 

where 

1 ( Pe _ a M^)2 
#i = 2 t P « + * 1 + 2R 2 . (50) 

/ D N(N-l) F„(Pu-aGu)\2 

2R ? r M 1"E M ^ 2 



(51) 

It is easy to see that for the special initial conditions 9^ = 0, H 2 = 0. Also the Poisson 
bracket {Hi,H 2 } oc H 2 vanishes for these initial conditions. The corresponding quantum 
statement would [H 1 ,H 2 \ oc H 2 , and therefore if we consider the subspace {tp} on which 
H 2 ip = then Hi will act invariantly on this subspace. Further H = Hi on this subspace 
and thus the eigenstates of Hi will be exact eigenstates of the full system and conversly. But 
the problem of solving Hi is analogous to that of two anyon problem with a — > Q , Ar ( i ^~ 1 ) _ 
Therefore these solutions are given by, 

N(N - 1) 

E = 2m+\j-a^^ J -\ + (N-l), (52) 

^ = CR lj - al exp- R2/2 f(R), (53) 
Hii) = Eip (54) 

where /(it!) is some polynomial of degree 2m in R. In general the normalisation constant 
C in ip will have dependence on 8^s, restricted by H 2 ip = 0. These are precisely all the 
known solutions and the above arguement shows that there are no more solutions satisfying 
Hiip = Eip and H 2 i\) = 0. Thus H 2 ip = characterises the subspace spanned by all the 
known exact solutions. So the quantum counterpart of classical partial separability implies 
the existence of these solutions. 
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V. INTEGRABILITY 



We now come back to the integrability aspect. In Sect. 3 we exhibited a set of 2(iV — 1) 
constants of motion and regarded the system as potentially integrable. However for the 



system to be integrable via action angle variables further conditions have to be satisfied [H 
If M(S a , I a) denote the set of points in phase space at which the constants of motion have 
values £ a , l a then (a) M is a 2(N — 1) dimensional submanifold iff £ a , l a are all independent, 
(b) if M is a submanifold which is compact and connected then M is a 2(N — 1) dimensional 
torus. If these conditions are satisfied then one can introduce action-angle variables in the 
neighbourhood of M. The orbits on M will all be conditionally periodic in general and for 
the periodic orbits one may use Bohr-Sommerfeld quantisation to get a subset of eigenvalues. 

We also saw in the previous section that there exist M(£ a ,Z a )'s for which action-angle 
variables can be introduced. These correspond to the collective motion of N-anyons and 
their semiclassical quantisation will reproduce the exactly known eigenvalues described ear- 
lier. The Hamiltonian in Eq.(50) being identical in structure to the oscillator Hamiltonian, 
the orbits of this Hamiltonian will be in one-to-one correspondence with the orbits of the 
oscillator and hence they will all be periodic. However if we go away from these initial 
conditions some of the orbits of oscillator will cease to be periodic when a dependence is 
introduced. For instance the Euler-Lagrange equations 



ij 



-fy (independent of a). 



For any pair ij, i ^ j, in general describes an ellipse in the configuration space. However 
for orbital angular momentum zero, the ellipse degenerates to a straight line, hence = 
is on such an orbit. Since these points are not admissible when a / these orbits cannot 
lead to periodic orbits in the phase space. For all choices of £ a , l a such that for some 
initial conditions the phase space orbits will have ry approaching arbitrarily close to zero, 
M(£ a , l a ) cannot be a torus. (M, even if a manifold will fail to be compact and/or connected.) 
Therefore for a ^ although we have a potentially integrable system it is not generically 
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integrable via action-angle variables. The existence of such an M indicates the possibility 
of "extreme sensitivity" to the initial conditions. This is a possible reason for the "level 
repulsions" seen numerically. This bears a resemblence to the billiard system considered by 
Richens and Berry [I3[] though the reasons for the failure of integrability via action-angle 
coordinates appears to be different. Following Richens and Berry we conclude that our 



system is Pseudointegrable 13 



VI. SUMMARY AND CONCLUSIONS 

In this paper we have exhibited two properties of the many anyon system: (1) The partial 
separability of collective motion and (2) pseudointegrability. We have shown that property 
(1) explains the existence of the exactly known eigenvalues which is some what uncommon 
for a generic many body problem. This also shows that the exactly known spectrum incor- 
porates only a some what trivial aspect of anyon dynamics. The nontrivial aspects although 
partially uncovered by numerical results are still elusive. Property (2) enables one to un- 
derstand qualitatively the origin of the "level repulsion" seen in some numerical results. 
Further study of pseudointegrability , for example the topology of level sets M, may enable 
one to adapt different techniques to get a handle on the nontrivial aspect of the spectrum. 
We have demonstrated that the N-anyon system is a nontrivial example of a many body 
pseudointegrable system apart from the known problem of a billiard in polygonal enclo- 
sures with a regular obstacle inside. In the light of the pseudointegrability, it seems that 
the characterisation of quantum integrability must be understood with further conditions , 
for example self-adjointness of the constants of motion in involutions and the existence of 
common dense domains. 

We thank R.Ramaswamy for drawing our attention to the notion of pseudointegrable 
systems and Radha Balakrishnan for discussions. 
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